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Abstract 

We study the local scaling properties associated with straight line 
periodic orbits in homogeneous Hamiltonian systems, whose stability 
undergoes repeated oscillations as a function of one parameter. We 
give strong evidence of local scaling of the Poincare section with expo- 
nents depending simply on the degree of homogeneity of the potential. 
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It has been recognized for some time now that periodic orbits 
play a crucial role, whether in classical or semiclassical dynamics. Of the 
infinite periodic orbits the most important ones are those with the shortest 
periods and highest stabilities. Gutzwiller's trace formula in the semiclassical 
quantization of chaotic systems [1], as well as the zeta function approaches 
in classical and semiclassical mechanics of such systems [2], accord them 
the highest weights. The bifurcation properties of these orbits also assume 
considerable significance. Certain atomic experiments have revealed the im- 
portance of bifurcations of closed orbits even when the dynamics is chaotic 
[3]. 

Most studies of classical Hamiltonian systems have focused upon single 
parameter systems, upon whose variation the system smoothly undergoes a 
transition from regular motion to chaotic motion via stages of mixed phase 
spaces. It can so happen that integrability may be suddenly recovered for 
certain values of the parameter. However, even while the parameter variation 
is over a range in which the remnant tori are being destroyed and replaced by 
chaotic trajectories, there could be periodic orbits that are rapidly undergo- 
ing stability oscillations implying the creation of secondary tori and regular 
regions in the phase space. This has been known for sometime and seems to 
be more generic with homogeneous Hamiltonian systems [4] . 

Such stability oscillations occur in very simple periodic orbits and as 
stated above these are of importance. Homogeneous Hamiltonian systems, 
while rather special, allow certain simplifications that make their study use- 
ful. While in general Hamiltonian systems the orbits form one parameter 
families with energy being the parameter [5], in homogeneous systems vary- 
ing of energy simply scales the orbits without changing the orbit structure in 
the phase space, that is bifurcations and related phenomena cannot occur as 
a function of energy, in general. Thus we resort to changing the Hamiltonian 
itself in the form of parameter variations. 

For homogeneous Hamiltonian systems, Yoshida [4] has given an exact 
and simple expression for the trace of the monodromy matrix of certain 
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straight line periodic orbits which have in general low periods and high sta- 
bility, or are among the least unstable orbits. The monodromy matrix is 
the linearization of the Poincare map in the neighborhood of the periodic 
orbit. To fix ideas and introduce the scaling associated with these orbits, we 
will begin with the well studied model of the quartic oscillator given by the 
Hamiltonian [6] 



where (3i, (32, a > 0. For fixed values of the parameters (3i^2, as a is increased, 
the phase space is known to become more chaotic. Although even at a = cx) 
there are islands of stability, these are of miniscule proportions. The straight 
line periodic orbit, the "channel orbit", specified by the initial conditions 
{Pi, Qi, 0, 0), is clearly one of the simplest orbits of the system and is known 
to play a crucial part in the semiclassics and quantum mechanics of the 
oscillator; for instance they scar a series of eigenfunctions which form a near 
WKB series even in the highly chaotic regimes [7]. 

The channel periodic orbits do not increasingly become unstable as a is 
increased, they recover stabihty by repeated oscillations. This imphes that 
over whole ranges of a, however chaotic the rest of phase space may be, 
there are islands of stability around this periodic orbit and that in the stable 
regions various bifurcations give rise to new periodic orbits. For the instance 
of the Hamiltonian specified by Eq. 1, the Yoshida formula [4] gives 



where J{a) is the monodromy matrix for the half Poincare map [8] of the 
oscillator. Thus J{a) is the linearized map about the channel periodic orbits 
specified by {q2 = P2 = 0). Due to the symmetries of the system we are 
considering, namely refiection symmetries about the various axes, the half 
map defined as successive intersections of the trajectories with the plane 
Qi = 0, irrespective of whether pi is positive or negative, is an one to one area 
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preserving map. We thus note that the orbit can change stabihty whenever 
a = f3im{l + 2m) where m is any integer, as at these values Tr J{a) = ±2. 

For large enough a, the phase space is mostly chaotic, hence when the 
channel orbit is stable, its island of stability must be rapidly shrinking with 
a. We can compare the stable areas on the half Poincare sections at various 
a, such that the central orbit stability is the same at these values, and the 
slope of the stability curve, dTr {J{a))/da, has the same sign. For instance 
Fig. 1 shows the neighborhood of the origin, corresponding to the channel 
orbit (g2 = P2 = 0), when its stability is just about to be lost in a pitchfork 
bifurcation, i.e. when Tr J{a) — 2 and the trace is increasing. It is clear 
that while the islands are shrinking with the parameter, they are essentially 
similar and would possibly scale with a. We thus formulate our principal 
results, which are at present only in the form of numerical explorations, as 
follows. 

Let 

(I2 = /(?2,P2; a) P2 = g{<l2,P2] ot) (3) 

be the half Poincare map. As a consequence of the reflection symmetries in 
the oscillator, the function / is such that f{q2,P2]o:) = — /(— 52, —P2', ct), with 
a similar relation for g. Then the scaling of the section implies the scahng 
of the above functions. Let us choose two values of the parameter a and a', 
such that say a' > a. If a and a' are related by Tr J(q!) = Tr J(q;'), and the 
stability is either increasing at both a and a', or decreasing, then 

(-r''f{q2,P2;a) = /((^)-7lg2,(^)-72^,;c'). (4) 

a a a 

Here 71 and 72 are the scaling exponents for the q2 and p2 directions respec- 
tively. 

For the class of Hamiltonians (we call here class I), given by 

U — "'",^2 I 1 2 I ^ {a „2n j_ o 2n\ . <^ / 2 2n-2 , 2 2n-2\ 
-H2n - -^Pl + 2^2 + ^(^^1^1 + -^252 ) + ^Wl92 + ^2^1 ): 
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of which the quartic oscillator used above is a special case, we conjecture the 
following, based on numerical evidence to be presented below: 

2n+l 2n-l 
^'^^ 

Thus for the Hamiltonian of Eq. 1, 71 = 5/8 and 72 = 3/8. One of the 
consequences of the above is that the area of the sections scale simply as 
a~^, independently of the degree of homogeneity of the potential. A similar 
scaling relation is found to be true for the function g{q2,P2] «)• 

The vahdity of the above scahng relationships are restricted to a certain 
region around the periodic orbit, in this case around the origin of the section. 
The scaling is in this sense only local. We have observed that the area 
of stability may be safely taken as the region in which the scaling holds, 
although this can be a serious underestimation, as will be shown below. 
We will illustrate the vahdity of the scahng by taking one of the outer most 
points of the stable region of the sections when Tr J{a) — 2 and is increasing. 
In this case there is one island chain consisting of eight islands, that have 
been earlier created, and have grown out and are near the chaotic sea (Fig. 
1). We will take the distances between the origin and the central period 
eight orbit to verify scaling. Let the period eight orbit's intersection with 
the positive 52 axes be at ^1(0;) and with the positive p2 axes be at d2{a). 
Fig. 2 shows the scaling of these distances with a, i.e., ^1(0;) ~ and 
d2{ct) ~ a~^'^. The lines shown are those of best fit. Their slopes are equal 
to —0.621 and —0.372 and are very close to those predicted by the above 
Eq. (5). The scaling seems to become better with increasing a, so that the 
first few points were neglected while calculating the slope. Increasing a leads 
to a deterioration of the accuracy of the numerical integrations. Hence we 
have used smaller step sizes of the order of 10^^ in a fifth order Runge-Kutta 
integrator for converging the exponents at these high parameter ranges. 

The exponents found from the above can be used to directly verify the 
scaling of the half first return maps as given by the Eq. 4. In the case of the 
Hamiltonian of Eq. 1 , Fig. 3 shows the absolute value of the difference of the 
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two sides of Eq. 4 for the function /, for the case when a' = 120 and a = 66. 
At these values of a the trace is 2.0 and increasing when the stabihty is about 
to be lost in a pitchfork bifurcation and there are large stable islands (Fig. 
1). Comparing figure 1(a) and 1(b) with Fig. 3 indicates that the area over 
which the scaling remains valid is much larger than the "area of stability" . 
A similar result is obtained in the case of the function g as well. 

To emphasize this we may take the case when Tr J{a) — 2 and decreasing, 
when the channel orbit is about to gain stability and create two new unstable 
orbits (for instance a' — 136 and a — 78) . In this case there is no stable 
island, yet the scaling of the functions is valid over a large range and the 
picture obtained is very close to that of Fig. 3. The scaling relation is found 
to be true even in the case when the central channel orbit is unstable. At 
this stage we note that the scaling of the functions / and g do not necessarily 
imply scaling of the orbits, as in a chaotic flow which is ergodic the phase 
points will explore regions in which the scaling is invalid. However, in the 
case when the orbits never leave the region of valid scaling, we can expect 
the orbits themselves to be scahng. This would happen if the central orbit 
were to be stable, and explains our interest in this range of parameter values, 
as well as the likeness in the sections of Fig. 1. 

Verifying scaling of the functions is much easier than measuring distances 
implied in Fig. 2. Using the exponents found when Tr J[a) — 2, we have 
verified using Eq. (4), the scaling laws with identical exponents independent 
of the value of the trace. Another rather efficient method of determining the 
exponents based on Eq. (4) is to assume a fixed initial condition with P2 = 
and searching along a range in the exponents for 71, as this is unaffected by 
the value of 72, and then searching for 72 using the 71 obtained from such a 
procedure. Identical scaling behaviour with the exponents given by Eq. (5) 
is observed when f3i 7^/32, i.e., when the symmetry of the above examples 
is broken into C2V, and this is illustrated also in Fig. 2 for the case f3i = 0.5 
and = 1.0, when the lines of best fit have slopes equal to -.622 and -.372. 

An almost identical picture is obtained when we take other class I sys- 
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tems. For instance we consider the Hamiltonians and Hs whose potential 
energies correspond to n = 3 and n = 4 respectively, within class I. The 
figures analogous to Fig. 2 is shown in Fig. 4 for these oscillators. The lines 
are once more those of best fit, and the slopes for the sextic are —0.589 and 
—0.422, while for the octic potential they are —0.566 and —0.442, which are 
very close to the values given by Eq. (5), when we consider that by definition 
7i,2 are negative of the slope. Once more we find that the scaling gets to be 
nearly perfect for large values of a. 

Potentials that contain terms which do not affect the stabihty of the 
channel orbits form different classes of Hamiltonians from those considered 
above. For instance, one simple set of Hamiltonians we call class H is of the 
form 

HL = \pI + \pI + ^(/^i??" + (^^^?) + ^«<?2 , n>2. 

The channel orbit is always marginally stable (Tr J {a) = 2), independent 
of a and there is a stable region around this orbit which continuously scales 
with a. We found the corresponding exponents to be well predicted by the 
following rule: 



so that the area still scales as only in the case when n = 4. The term 
qiq2 {n > 2) is like a "gauge term" as far as the stability of the central orbit 
is concerned. In this class of Hamiltonians the symmetry of parity is broken 
when n is odd, and the potential in these cases is bounded only if — 1 < o; < 1. 
If n is odd, the half Poincare map defined earlier for class I Hamiltonians is 
not valid, and hence we use the usual definition of full Poincare map, namely, 
as the successive intersections of the trajectory with the plane gi = and 
Pi > 0. For example, in case of the Hamiltonian specified by the potential 
{qf + q2)/Q + the phase space is largely chaotic for the seemingly 

low values of the coupling parameters near unity. Complete chaos is however 
absent, not only because of the channel orbit but also due to the existence 
of one more stable island. In this case the exponents were found, using the 
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methods specified above, to be 71 = 72 = 1. The generahzation, Eq. 6, is 
based on similar computations for larger values of n (up to n = 7). 

We have briefly noted some, what we believe are new, local scaling be- 
haviours of certain homogeneous Hamiltonian systems. The above being 
in the nature of preliminary numerical exploration, we cannot exhaustively 
comment on the classes of Hamiltonian systems with distinct scaling laws, 
even within the sub-class of homogeneous systems. The number of degrees of 
freedom we have considered in this Letter is only two and higher dimensional 
generalizations while interesting have not yet been explored. It is also not 
clear if such scaling behaviours can be observed in non-homogeneous sys- 
tems with similar periodic orbits. In future work we hope to address some 
of these questions as well as study the semiclassical implications, if any, of 
such scaling. 
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Figure Captions 



Figure 1 Poincare surface of section around the origin for the quartic oscillator 
cases a) q; 66 and h) a = 120, with Pi — P2 — 

Figure 2 Scaling of the distances for the quartic oscillator, when a' — 120 and 
a — 66, and a similar case when j3i = .5 is also shown. The upper two 
and lower two lines correspond to ^2(0;) and di{a) respectively. 

Figure 3 The absolute value of the difference between the left and right hand 
sides of Eq. 4 for the case when Tr J{a) — 2 and increasing, a' — 120 
and a — 66. 

Figure 4 Scaling of the distances for the case of the a) sextic and b) octic oscil- 
lators (/3i,2 = 1.0). 
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